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ABSTRACT 

We have obtained the mass-metallicity (M-Z) relation at different lookback times for 
the same set of galaxies from the Sloan Digital Sky Survey, using the stellar metal- 
licities estimated with our spectral synthesis code starlight. We have found that 
this relation steepens and spans a wider range in both mass and metallicity at higher 
redshifts. We have modeled the time evolution of stellar metallicity with a closed- 
box chemical evolution model, for galaxies of different types and masses. Our results 
suggest that the M~Z relation for galaxies with present-day stellar masses down to 
10^°Mq is mainly dr iven by the history of star formation history and not by inflows 
or outflows. 
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1 INTRODUCTION 

Since the study by iLegueux et al.l l| 19791 ). who found 
a luminosity-metallicity relation for irregular galaxies, 
many papers have reported the existence of a luminosity- 
metallicity or mass-metallicity {M-Z) relation for all kinds 
of galaxies: The more massive galaxies are also the ones with 
more metal-rich interstellar medium (ISM). 

Considerable work has also been devoted to explaining 
the M-Z relation. For instance, Trcmonti ct al. (2004) find 
that outflows play an important role in shaping the M-Z 
relation, as the weaker gravitation potential well of less mas- 
sive galaxies makes them more prone to los e enriched gas via 
galact ic winds or supernovae explosions. iFinlator fc Davj 
l|2008l ) claim that inflows of pristine gas could also be a 
sound explanation, as the same amount of gas falling into 
galaxies would have a greate r impact on the chem ical com- 
position of less massive ones. iKoppen et all l|2007l ) advocate 
the role of an integrated initial mass function: Massive galax- 
ies comprise more massive clusters which may contain more 
meissive stars that enrich the ISM faster. 

The M-Z relation derived fro m emission lines has 
been found to chan ge with redshift (jSavaglio et al ] I2OO5I : 
iMaiolino et al . 2008). The comparison between M-Z at low 
and high redshifts is, however, not straightforward. Differ- 
ent populations of galaxies live at different redshifts, and 
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this alone can affect the derivation of the nebular metallic- 
ity, even if one were to apply the same calibrations for both 
samples (seeiStasihska 2008). Nevertheless, at least qualita- 
tively, the observed results are in agreement with what is 
expected for the evolution of galaxies. 

This brings us to another concern in M-Z studies: the 
bias in the derivation of the nebular metallicity from emis- 
sion lines. The evolution of the mass-metallicity relation is 
usually studied in terms of th e abundance of oxygen i n the 
ISM gas. As known (see e.g. iKewlev fc EllisonI [20oi ) . the 
measurement of nebular abundances is very dependent on 
the method and calibrations used. 

With our stellar population synthesis code starlight, 
we can determine the stellar metallicities, Z^, and the to- 
tal mas ses in stars, M-^, at differ ent epochs for any given 
galaxy (|Cid Fernandes et al.|[2007l ). We can thus follow the 
evolution of the M-Z re lation for the same set of galax- 
ies at different redshifts. iPanter et al.l (|2008| ) used similar 
techniques and discussed the nebular metallicities with re- 
spect to stellar metallicities. In this study, we rely on stellar 
metallicities alone because, although they have their own 
problems, they are free from the biases affecting nebular 
abundance determinations. The use of Z^ also allows us to 
explore more massive galaxies, for which the nebular metal- 
licity estimates are not possible, such as AGN hosts and 
galaxies which have stopped forming stars. 

Since starlight allows the determination of the star 
formation histories (SFHs), we can go one step further and 
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Figure 1. Histories of our parent sample in stellar mass bins. 
Thicker lines are used for more massive bins. Panels show the 
median curve for time evolution of stellar mass evolution (left; 
see Eq. [TJ and mctallicity evolution (right). 



construct a simple chemical evolution model using these 
SFHs and see under what co nditions it reproduce s the ob- 
served time evolution of Z.. ISavaglio et all l|2005h made a 
similar study, but using a theoretical SFH with e-folding 
time, related to galaxy mass, and for nebular metallicity in- 
stead of . Here we use the SFH obtained directly for each 
galaxy. 



2 DATA 

Our parent sample is composed of all objects spectroscop- 
ically classified as galaxies from S*^** Data Release of the 
Sloan Digital Sky Survey (SPSS') (|Adelman-McCarthv et all 
I2OO7I ). We exclude duplicate observations, and impose the 
following selection criteria: 14.5 < rrir < 17.77 (from the 
Main Galaxy Sample definition), z-band aperture covering 
factor > 20% (to reduce aperture effects and avoid intra- 
galactic sources), S/N > 10 at 4750 A (to provide reason- 
able stellar population fits), and a narrow range in redshift 
2 — 0.11 < 0.015 (to be able to transform stellar ages into 
redshift in a simple way). The central redshift of 0.1 was 
chosen to maximize the number of objects and thus ensure 
reliable statistics over the whole range of ages considered in 
the analysis. We adopt a Ho = 70kms~^ Mpc"'^, Q,m = 0.3 
and Qa = 0.7 cosmology where appropriate. There are 82662 
objects in the resulting sample. 

We use the same d ata processing as outlin ed in 
ICid Fernandes et al.1 (|2005l ) and iMateus et all l|2006l ). Our 
code STARLIGHT recovers the stellar population content of a 
galaxy by fitting a pixel-by-pixel model to the spectral con- 
tinuum. This model is a linear combin ation of 150 simple 
stellar populations (SSP) extracted from lBruzual fc CharlotI 
l|2003l . hereafter BC03) with ages 1 My r < U < 18 Gyr, and 
metal licities 0.005 < Z/Zq < 2.5, as in lCid Fernandes et al.l 
Emission lines are then measured in the residual 
spectra, which reduces the contamination by stellar absorp- 
tion features. 

We divide our parent sample into three sub-samples: 
star- forming (SF), retired (R) and Seyfert (S) galax- 
ies. For both the SF and S samples, we select objects 
with S/N > 3 in the [O iii]A5007, H/3, [N ii]A6584 
and Hq emission lines, i.e., a l l the lines involved in the 
[Baldwin. Phillips, fc TerlevichI (|l98ll . BPT) diagram. The 
SF sample is composed by the objects below the line pro- 



posed bv lStasihska et al.l (|2006l ) to separate SF galaxies from 
AGN hosts. The Seyfert sample is defined as the galax- 
ies with lo g([N ll]/Ha ) > —0.3 and above the line de- 
fined by iSchlickmannI (|2008l ) in the BPT diagram to sep- 
arate Seyferts from LINER-like galaxies: log([0 lll]/H/3) = 
0.90 log([N ii]/Ha) 0.48. 

The R sample is composed of galaxies with little evi- 
dence of either nuclear or star-forming activity, which in- 
cludes passive galaxies, galaxies with very weak emission 
lines and even some which would be trad itionally classi- 
fied as LINERs (see IStasiriska et al.1 lioO^) . Therefore, we 
impose the following cuts for the R sample: S/N < 3 in 
the 4 BPT lines and [O n]A3727 ("passive"); or S/N < 3 
in [O III] and/or H/3, but S/N > 3 in [N 11] and Ha, and 
log([Nll]/HQ) > -0.15 ("weak emission lines"); or S/N > 3 
in the 4 B PT lines, log( [N ii]/Hq) > -0.15, below the line 
defined by ISchlickmann ! (2008), and Ha emission predicted 
by old stellar populations equal or greater than the observed 
value ("retired mimicking LINERs") 

We have also divided our sample and sub-samples into 
six present-day stellar mass bins centered in log M^/M© = 
10.0 (A), 10.3 (B), 10.6 (C), 10.9 (D), 11.2 (E) and 11.5 
(F), each on e 0.30 dex wide. In the same vein as our pre - 
vious works (|Cid Fernandes et al ][2003; I Asari et al.| [2007ll. 
defining bins allows us to alleviate the problem of the mass 
completeness of the sample. We consider narrow mass bins 
in order to have very similar galaxies inside each bin. Al- 
though the starting sample contained a small proportion of 
galaxies with masses smaller than 10^° M©, we do not con- 
sider them in the current analysis. Indeed, due to the redshift 
cut we imposed, only the most luminous low-mass objects 
are present and they are not representative of the bulk of 
galaxies of their mass range. 



3 GALAXY HISTORIES IN MASS BINS 



As explained in ICid Fernandes et al.l (|2007l '). starlight re- 
covers the fraction Xj that each SSP j contributes to the 
total fight of a galaxy in the spectral range covered by the 
SDSS. We can translate Xj into the mass fraction presently 
locked into stars, fij , or the fraction of mass ever converted 
into stars, fij. 

The stellar mass history can be followed as a function 
of lookback time t by summing up the converted-mass frac- 
tions: 



(1) 



The mean stellar metallicity at a lookback time t is 
defined as the total mass in metals locked in stars divided 
by the total stellar mass at a given time: 



z4t) = Y.fijZ, 



(2) 



Another way to look at the metallicity evolution would 
be to consider the average metallicity of all the stars born 
at the same epoch, for each value of t. This, however, is a 
much less robust quantity, and we found that SSP bases at 
hand nowadays do not allow such detailed definition. 

Figure [1] shows the median values of ry* and versus t 
for mass bins A-F of our sample. As previous studies have 
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Figure 2. Left: Present-day mass-metallicity relation for all the galaxies in our sample. The thick line marks the median and thin lines 
the quartiles. Right: Mass-metallicity evolution for each of our present-day stellar mass bins A-F (black). From top to bottom, bold lines 
mark the mass-metallicity relation for redshifts z = 0.1 (t = Gyr), 0.2 (1 Gyr), 0.7 (5 Gyr) and 1.9 (9 Gyr). 



alrea dy shown (e.g. ICid Fernandas et al ] |2007l : IPante7 et al.l 
I2OO8I . for SF galaxies only - here we extend the study for all 
galaxies), there is a clear downsizing effect: The less massive 
galaxies are slower in converting their mass into stars and 
in building up their metal content. 



4 THE OBSERVED M-Z EVOLUTION 

Figure [2] shows the mass-metallicity relation. The left panel 
shows the present-day M-Z relation for all the galaxies in 
our sample. The right panel depicts the evolution of the M- 
Z relation for the galaxies in mass bins A-F (black lines): 
Zi,{t) vs. Mi,{t), where Mi,{t) is the stellar mass integrated 
from the epoch of galaxy formation until a lookback time t, 
assuming that the stars always pertain to the galaxies where 
they were born. The resulting M~Z relations for different 
redshifts are plotted in bold lines. One can see that the mass- 
metallicity relation evolves. 

As lookback time increases, the M-Z relation steep- 
ens and cove rs a larger range of values. As investigated by 
iPanter et al] (t2008). using stellar modelling codes other than 
that of BC03 changes the values of M*(i) and Zi,{t), but 
the general behavio ur seen in the present paper as well as in 
iPanter et all l|2008l ) remains. One should note that the M- 
Z relation and evolution is significantly sample-dependent, 
which makes difficult any quantitative comparison between 
works by different authors. 



5 A SIMPLE CHEMICAL EVOLUTION 
MODEL 

We now investigate the behaviour of the M*-Z* relations 
shown in Fig. [2] with the help of a simple model of chemical 
evolution of galaxies. We assume a closed box with initial 
metallicity equal to zero, and use an instantaneous recycling 
approximation and a constant yield. The mass-weighted av- 
erage mass fraction of metals at lookback time t is then given 



by (Edmunds 1990, eq. 35): 
/(f)ln(l//(t)) 



where f{t) is the gas mass fraction and y is the ratio of the 
mass of metals released per stellar generation to the total 
stellar mass locked in remnants. The gas mass fraction is 
given by: 



fit) 



1 - R 
Mt 



SFR(t') dt', 



(4) 



where Mt is the total mass of the galaxy (stars and gas), 
SFR is the star formation rate, R is the returned mass frac- 
tion due to stellar winds and supernovae and to is the age of 
the oldest SSP in the base. The star formation rate SFR(t) 
can be obtained from the results of the starlight synthesis: 

S¥R{t)dt = M^^il{t), (5) 

where is the total mass converted into stars (i\f^ = 
Mi, /{I — R)). Using Eqs. [1] and [5l f{t) becomes simply: 



fit) 



(1 - fno^)v*it), 



(6) 



Z4t) = y^y- 



(1 - fit)) 



(3) 



where /now is the present-day gas mass fraction. Knowing 
the value of y, one can use Eqs. |3] and [S] with rji,{t) given by 
STARLIGHT to find what value of /now is needed to reproduce 
the observed evolution of the mean stellar metallicity Zi,{t). 

Unfortunately, y is not well known. The measured stel- 
lar metallicities are mainly sensitive to the abundance of 
iron, and the Fe/O ratio is expected to vary with time. For 
our simple approach, we assume that all the metals vary 
in lockstep and simply speak of the "metal yield". We de- 
termine an empirical value of y by applying our model to 
the highest ma ss bin of the SF sample, using the Schmidt- 
Kennicutt law (|Kennicuttlll998l 'FI to derive the value of /now 
from the extinction-corrected Ha luminosity. This procedure 
should give a reasonable value of the metal yield, since it is 
for the highest mass bin that the closed box model is ex- 
pected to be the most relevant. We obtain y — 0.03. This 
value is in satisfactory agreement with the value (0.04) ob- 
tained for a iChabrier. (,2003 ) stellar initial mass function by 



^ We use the Schmidt-Kennicutt law for nor mal disk galaxies 
(biva riate least - squar es regression in fig. 2 of iKennicutg Il998t) 
for a IChabrjd |20p3) initial mass function, and including the 
contribution of helium in the gas mass fraction. 
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Figure 3. Chemical evolution in our samples of SF, R and S galaxies in the 0.30 dex-wide stellar mass bins, centered in log M^/Mq = 10.0 
(A), 10.3 (B), 10.6 (C), 10.9 (D), 11.2 (E) and 11.5 (F). Each panel shows the median and quartiles of the evolution of as found by 
STARLIGHT (black lines) and as obtained with the simple closed-box model (red lines). The value of /now needed to reproduce the median 
present-day Zi, is indicated at the top right of each panel. For the SF sample, we also indicate the median of /now as measured by the 
Schmidt-Kennicutt law (/sk). The yield y was adjusted to reproduce bin E of the SF sample. The number of galaxies in each panel is 
indicated at the bottom. Panels were left empty if they contained less than 250 galaxies. 



multiplying the oxygen yield from lWooslev fc Weaveij (|l995l ) 
by the Solar metal-to-oxygen mass fraction. 

We therefore adopt y = 0.03 and apply our model to 
try to reproduce the observed metallicity evolution for all 
the mass bins. Figure [3] shows Zi, (t) as a function of t for 
bins A-F of our SF, R and S galaxy samples. The black lines 
represent the median and quartiles of the observed distribu- 
tions. The red lines are the model results for the values of 
/now indicated at the top of each panel. The models were 
adjusted by eye to reproduce well at least the last 1 Gyr of 
the observed Zi,{t). For the SF sample, we also indicate the 
median value of /now obtained from the Schmidt-Kennicutt 
law. This figure shows that this simple model can reproduce 
the observed metallicity evolution reasonably well, given the 
crudeness of our approach. The discrepancies for ages above 
10^ yr are at least partly due to t he known problems faced by 
spect ral synthesis in this range (|Gomesll2005l :l Koleva et al.l 
|2008| ). In particular, the jumps in both rif, and at ~ 1 
Gyr ar e artifacts related to the spectral libraries used in 
BC03 I Cid Fernandes et all |2008| ). New evolutionary syn- 
thesis models with more accurate stellar spectra should im- 
prove the recovery of SFHs in this range. 

It is worth of noting that /now is at least qualitatively 
in agreement with expectations: it is larger for SF galax- 
ies than for retired galaxies, and decreases with increasing 
M*. As a matter of fact, for all the mass bins of the SF 
sample /now is compatible with the values derived from the 
Schmidt-Kennicutt law (bin E is compatible by construc- 
tion) . Seyferts also shows little residual gas mass. The chem- 
ical evolution of Seyferts is (within our resolution) similar 
to the evolution of non-Seyfert galaxies. 



To sum up, for present-day stellar masses down to 
10^" M0 and within the current limitations of astropaleon- 
tology, the closed box model explains the present-day M-Z 
relation and its evolution quite well if one uses the SFR(t) 
obtained directly from the spectral synthesis of each galaxy. 
This is not to say that infall or outflow do not play a role. 
What we want to stress is that it is the star formation his- 
tory which is the main driver of the mass-metallicity relation 
and its evolution. 



6 SUMMARY AND DISCUSSION 

Until recently, the M-Z relation in galaxies had been mostly 
studied using abundances derived from emission lines. Ex- 
ploiting the stellar metallicities obtained with our spectral 
synthesis code starlight, we have obtained the M-Z rela- 
tion at different lookback times for the same set of galaxies. 
The use of stellar metallicities, even if those are not very 
accurate, has several merits. It allows one to probe metal- 
licities at different epochs of a galaxy evolution. It makes it 
possible to study galaxies in a larger mass range. Galaxies 
without emission lines or galaxies with an active nucleus do 
not have to be discarded. Finally, with stellar metallicities, 
one avoids the biases that affect the statistical methods used 
to derive metallicities from emission lines. 

Our main results are the following. We have found 
that the M-Z relation steepens and spans a wider range 
in both mass and metallicity at higher redshifts. The more 
massive galaxies show very little evolution since a look- 
back time of 9 Gyr. This means that they evolved very 
fast in a distant past, supporting the downsizing scenario. 
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This is in agreement with other studies of the build-up of 
metals in SDSS galaxies a s rev ealed by fossil recor ds, e.g., 
ICid Fernandes et all l|200'it and lPanter et al.l ([2008'). 

We have modeled the observed time evolution of the 
mean stellar metallicity using a closed-box, instantaneous 
recycling chemical evolution model, for galaxies of different 
types and masses. We find that this model is compatible 
with the observations. This suggests that, in the mass range 
studied here (from log M^/Mq = 9.85 to 11.65), the M-Z 
relation is mainly driven by the history of star formation 
and not by inflows or outflows. By compa ring the nebu- 
lar an d stellar metallicities of SDSS galaxies, iGallazzi et al.l 
had argued that galaxies are not well approximated 
by closed-box models and that winds may be important. 
Our study suggests that those processes are not dominant, 
although they may play a role in the scatter of the M-Z 
relation in our mass range. 

Low-mass galaxies are not present in our sample be- 
cause of our selection in redshift. They would require a 
dedicated stud y to b e co mpared t o the result found by 
iTremonti et all l|2004 ) and lGarnetl l|2002ll that such galax- 
ies do suffer strong outflows. As a matter of fact, Garnett 
finds that outflows are important for galaxies with _B-band 
magnitude > —18. Only 6 out of the 82662 objects in our 
sample have g-band magnitude > —18. Using the Schmidt- 
Kennicutt law to compute the effective yields for each galaxy 
(i.e., the yield calculated from the obs erved nebular m e tallic - 
ities assuming a closed box model), ITremonti et al.l ()2004l l 
concluded that galactic winds are ubiquitous and very ef- 
fective in removing metals from galaxies. This however con- 
cerns only galaxies with masses up to IQ^^ Mq, which are not 
present in our sample. Our work emphasizes that there is a 
large range of galaxy masses - almost two decades - where 
the closed box model seems to give a valid representation of 
the chemical evolution of galaxies. 

Our models are parametrized by the present-day gas 
mass fraction, /now. We have found that the values of /now 
that allow us to reproduce the observed Zi,{t) histories are 
larger for less massive galaxies, in agreement with expec- 
tations. For galaxies of similar masses, they are larger for 
star-forming galaxies than for retired galaxies. The values 
of /now we flnd for SF galaxies are compatible with those 
obtained from the Schmidt-Kennicutt law. This sound re- 
sult supports the validity of our interpretation. 

Within our uncertainties, Seyfert galaxies have as little 
gas as retired galaxies. 
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